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DYNAMICAL THEORY OF STABILITY FOR ELASTIC RODS 
WITH NONLINEAR CURVATURE AND TWIST 

J. Wauer 

SUMMARY /l6l* 


Considering non-linear terms In the curvature as well as In 
the twist the governing boundary value problem for lateral ben- 
ding of elastic, transverse loaded rods are formulated by means 
of Hamilton’s principle. Using the method of small vibrations, 
the associated linearized equations of stability are derived, 
which complete the currently accepted relations. The example 
of the simplest lateral bending problem Illustrate the Improved 
effect of the proposed equations. 


1 . PROBLEM FORMULATION 

The stability equations which describe the tipping of elastic 
rods have been derived several times recently [1, 2, 3, with 
the assumption that the deviation between the undeformed and the 
deformed rod axis Is so small that differentiations with respect 
to the two Important position coordinates can be considered to 
be identical. 

This restriction will be given up here. The Hamilton principle 
described In [^] will be the point of departure for our analysis. 
Therefore, It Is the purpose of this paper to formulate this 
variational principle for describing the kinetic tipping with 
more extensive consideration of the rod curvature during defor- 
mation. Synthetic equilibrium analyses, which Is the continuation 

“Numbers In margin Indicate pagination In original foreign text. 
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of the Federhofer analysis, were performed early by Chwalla [5]. 
However, these have been forpiotten. It was only Ziegler [6] who 
considered this effect In the case of a compressed and twisted 
circular cross section. However, he did not pursue the most 
Important aspects of this process. In our paper. It will also 
not be possible to derive the consequences for the most general 
tipping problem of elastic rods with longitudinal and transverse 
loading. We have to restrict ourselves to cases where the rod 
axis does not experience any longitudinal extension. 

Using the Hamilton principle, we obtain the nonlinear boundary 
value problem from a variation. After this, one can transfer to 
the linearized stability equations using the method of small 
oscillations . 

In order to focus the reader's attention to the Important 
aspects of the problem, we will consider the Prandtl tipping 
rod with a pure moment load, without referring to a concrete 
special case. 


2. VARIATIONAL PRINCIPLE 

Let us consider an originally straight rod without preliminary /l62 
deformations of length 1. We assume that It has the constant mal*. 
bending stiffnesses EJ., EJ^ (E modulus of elasticity, 2 

moments of Inertia), the constant twisting stiffness OJ^ (0 shear 
modulus, J,p torsion area moment) and the constant mass per unit 
of length ((> Density, A cross-section area.) 

V/e will assume all of the usual assumptions made in classical 
technical bending theory, and assume pure Salnt-Venant torsion. 

We assume that the sheai* center and the center of gravity coincide. 
Also we will not consider the Influences of damping. 

The principle of Hamilton 

A / (•/• -Vi - »'.) (it + /AW dt^o ( 1 ) 

Is convenient! used for the systematic analysis. In order to 
evaluate it, we therefore must fl. st determine the kinetic energy 
T, the elastic potential , the potential of the external forces 
V^, and the virtual work 6V/ of the potential-free forces. 




( 2 ) 


The kinetic energy Is given by 

r- ;e/(//rj.rf^|rf*. 

where ff Is the absolute velocity of a material rod particle P; 

S, points along the arc length of the rod which has a curvature In 
Its general position. 

The elastic potential Is calculated from the sum of the 
shape change »nergy for bending 

y,B - { ajx fp \ + \ Ej, fp\ ( 3 ) 

o • 

and for twisting, 

O 

where p, - are the curvatures and q Is the torsion of the rod axis 
[7. 8], ’ 

These relationships can be generalized for the more general 
tipping process for both transverse and axial loading, by adding 
the following longitudinal extension energy 

(5) 

O 

within the elastic potential Vj , where e Is the longitudinal ex- 
tension of the central line of the rod. 

We will discuss the calculation of V„ and 6W later. 

a 


3. ANALYSIS OF THE STATE OF DEFORMATION 

Two coordinate systems shown In Figure 1 are used to describe 
the rod deformations. A first, space-fixed x, y, z reference 
system (unit vectors e, ) Is attached to one of the rod ends, 

where the z-axls coincides with the undeformed rod axis. A second 
body-fixed C, 1, C system (unit vectors ) Is established 

at an arbitrary point z. It points along the main moments of In- 
ertia of a cross-section and is rotated with respect to the space- 
fixed system, and It Is also tangential to the central line of the 
rod. 


Figure 1: Deformation state of a rod, cross section. 
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Figure 2: Description of’ deTormatlon using Euler angles. 

Consequently, the coordinate origin S' of the body-fixed system 
is therefore the geometric locus of all deformed cross-sectional 
centers of gravity. In the general case, It Is displaced by the 
coordinates u, v, with respect to the Initial point S. These 
quantities are measured In the space-fixed x, y, z system. The 
second coordinate system, whose rotation with respect to the 
(x, y, z) system. Is described using the Euler angles v'.O.tp ac- 
cording to figure 2, describes the stiffnesses EJ. « ard GJ . 

i ^ C 1 

Primes will refer to derivatives with respect to the arc length 
S. From the Klrchhoff gyroscope analogy [9], one finds the 
kinematic relationships 

R tp’ sin sin -4 cos ip, =» y' sin cos y — sin y, ^ = y' cos <p' ( 6 ) 

between the characteristics p^^ 9 I’od axis and the Euler 

angles v. 9 . 

As can be seen from Figure 2, the differential quotients 

*' •n sin sin y, y' = —sin t? cos y , z' *= cos ( 7 ) 

relate the space-fixed coordinates x, y, z and the Euler angles 

y. r . 

Tills means we are now In the position of formulating the 
Hamilton principle (1) as a function of tho de‘formatlon quanti- 
ties V- V" . In order to specify the boundary conditions, however, 
it is necessary to transfer from the Euler angles y.O.f to the 
center of gravity displacements u, v and the effective rotation 
of a cross-section t. 


We will assume small deformations, and therefore In (6) and 
(7) we will assume a small Inclination of the central line of 
the rod with respect to the space-fixed z-axls. However, we will 
not linearize with respect to v) as was done In [2], but we will 
expand up to quadratic terms. Then, from (6), we obtain 

v'o *in 4- co% 9 f«, p, wm y'fl %\tif , ^ y' •*- T* (8) 

and from (7) we find 

x'-Osinv, y»m^ftco%Y, I'-il— (9) 

The two first equations In (9) also apply for the derivatives of 
the dependent variables u, v In the x and y directions [10], If 
we differentiate with respect to the arc lenptth , we then find 

m" ■■ 0' ^in 4 co« y , V" —0' co» y 4- y . ( 10 ) 

If we formulate the expressions 

u” «ln (y +’ r) — v" (y 4- r) , co» (y 4- y) 4- »" «ln (y 4- y) ,) (11) 

(V4 y)' - ’ (.4V'_tV'). • .( 

we can see that they coincide with p^ 2 » 9 according to (6) If 
we use an addition theorem, where we take Into account the "Quad- 
ratic approximation". According to figure 2, the effective ro- 
tation angle of a rod cross-section angle Is 

T - y 4 y (12) 


for small Inclinations. Therefore, the curvatures p^ 2 
torsion q can be written In the form 

a u" sin T — n” rn% r , f t •* ^ t •**' **** ^ • (13) 

y-f- :(iiV'-rV') (li,) 


which only contain the usual quantities u, v, t and their deriva- 
tives . 

If we also assume a small twist angle t, then In addition to 
the torsion q (1^) we can write the curvatures p^ ^ In the fol- 
lowing form 

4 T«" . /k, - I." 4 TI-" (15) 


up to and Including quadratic coupling terms. Therefore, we have a 
consistent approximation for the quantities p^ 9 which are based 

on the theory of elasticity, and which Is not found In the literature 


5 


according to research by the author. There, one usually finds 

/*, . -p" » T«" . p, - m" 4 rp" . V - r' (16) 

[2, ^]. This Is a consequence of the Initially formulated assump- 
tion that derivatives with respect to a and z can be considered to 
be Identical. The coupling terms In (14), (15) have the same order 
of magnitude, and therefore only the approximation given here will 
be completely satisfactory, even though In practice the original 
relationship (l6) will remain useful. 

When calculating the absolute velocity, there are no special 
features compared with [4]. We find 


**p — (« — fir) e, -f (i» + (r) , 


(17) 


where C, n are the coordinates of P* In the local systems and dots 
are derivatives with respect to time. 

This means that the potential and the kinetic energy T for 
the tipping support can be formulated without the longitudinal 
extension of the rod axis, and without even considering a specific 
problem. If we Introduce the radius of Inertia k 

s 

I 

r - J (i* 4 - »• 4 - **t») * (18) 


which Is known from the literature [4], as well as 

j I • 

r, + tm”)** 4- v/ («" 4-rf-")*rf< 4--^- 


(19) 


If (14), (15) are substituted In (3), (4). 


4. TIPPING ROD WITH MOMENT LOAD 

In order to calculate the external potential In the virtual 
work 6W, one must specify the problem. We will consider a bending 
rod loaded by discrete end moments. We will not consider contin- 
uous bending moments, or continuously distributed forces and 
torsion moments. 

In addition to the constant direction moment pairs 

<=» o) « »' /) ” —Af.e, ( 2C ) 

we will assume the associated loads 
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• o) . A/^,(0 . 3f^(t - - .U,r,(i) 


( 21 ) 

/U)5 


Together they bend the rod In the bnslc state around the prin- 
cipal axis which Is stlffer (bending stiffness EJ^), before the 
stability limit Is reached. 

According to the Klrchhoff analogy, we car. then use the rela- 
t lonshlp 

till ••(.'/,/>, ♦ (22) 

for calculating the work [9]. The curvatures ^ torsion 

q are known as a function of u, v, t from (1^4), (15). M;., M , M 

are the components of the moment vectors (20), (21) in the f; , n, C 
system. Therefore, we must only decompose the constant direction 
moments (20) Into body-fixed coordinates. If we carry out this 
decomposition, then Mg can be written In the form 

,»f^(o) - M,(r, - rr, ♦ . M^{1) - -+ rr, - iiVj)!, ( 23 ) 

and according to (22), we finally get 

Air « Af,(Ar' - u' At)(1 ♦ M,{Av' - t 4M’)(i . ( 2^1 ) 

It Is found that a moment load which follows the rod as well as a 
constant direction moment load does not have a potential. Therefore, 
we have 

(25) 

and the virtual work Is finally given by (2^). 

The variations specified In (1), can then be evaluated using 
(l8), (19), (2^), and (25). The results are the coupled and non- 
linear motion equations for tipping In space of an elastic rod, 
with constant direction and following mom:?nt loads and nonlinear 
boundary conditions. We find the follow! ig for the bending defor- 
mations 

F . + {FJ, - EJ,) {rt.")” 4-^f(rV')' + + /i« - o . 

GJr ’ i 2 f>) 

EJ,v'- + (EU - EJ^) (TM ')" -^\{r’u"Y + + /ir - o 


where 


7 


Fjy' + (FJt - FJti rr ’ 4- + M,r - o 

or 

N'i» 0 • 

Fjy’ 4 (FJ, - £J,) (ny - o 

or 

H - 0, 

FJ,P" 4 (FJ, - T'w'- (M, M,) o 

or 

r'-o, 

Fjy" 4 (FJ, - FJ,) (tuy - (rwy - o 

or 

0 

1 


For the twist, we find 

- i - (^/i - F.U u’V - /!*> .0 < 2 8 ) 

where 

CtJj\r' — 1 (wV' — -4 My "O or T — o, ««o «nd / (29) 


Here we will Investigate the stability basic state which here 
Is a pure bending deformation around the stlffest bending 
principle axes. On top of this, we superimpose the Infinitesimal 
additional displacements 

MaaO-fil, r«»f, + 5, T — Oft (30) 


which are the components of the basic motion. These trial solu- 
tions (30) are substituted Into the nonlinear boundary value prob- 
lem ( 26 ) to ( 29 ). For the zero-order terms In the u, v, x we 
find a boundary value problem whose solution can be given In the 
form 


«« 

»’o 


M, + M , 


conM. , 


M, * M, / I 
iJx V “ * 




(31) 


The first order terms In the dashed quantities (we leave the dash 
off) results In the differential equations /166 


Kjy- + (A7i - F.J,) ^ + ,xu 

F’Jr [»' - T (»'>' - »>")]* - (Fh - F.h) r’y - fik*,r 


o.| 

It I 


and the boundary conditions 


FJ,W ■(- [(FJ, - FJ,) V- + M,] T 4- .Ct' - 0 

or 

m' « 0. 

FJ,h"‘ 4- (FJ, - FJ,) ^ (.•>')' - 0 

or 

W « 0 , 

F>Jt [x' - 7 (»'oM' - »C«")| f My - o 

or 

T ■■ 0 , 


(32) 


( 33 ) 


which are the linear stability equations. A comparison shows 
that the nonlinear Inl’ lal problem (26) - (29) ana the stability 
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equations (32)* (33) are a furtner continuation of the relationships 
Riven In [^), caused by the torsion q which has been expanded by 
the nonlinear terms q (1^4), 

As a specific example, we will first consider a rod which is 
supported by forks. Just like a Joint, such attachment cannot 
take up any bending moment j however. It does prevent the transverse 
displacement and the twisting of the end cross-section. The stabi- 
lity equations are substantially simplified and the Interesting 
difference between the following and the constant direction end 
moments Is lost. Both are conservative and we will simply write 
M for their sum (M^, + M^) In the future. For the theory, this 
deficiency Is not very serious, because the more complete consi- 
deration of rod torsion (1^) Is still In effect. There only exist 
a number of experimental r'' suits for a tipping support supported 
by forks. 

This means that the Inertia terms In (32), (33) do not have to 
be carried alonfe. for this special case. The tipping moment of 
Interest is the smallest load In which the non-trlvlal equilibrium 
position exl.'^ts. This means that It Is the lowest branching point 
of the eigen value problem (32), (33) which does not Involve time. 

Because of the fact that the coefficients partially depend on 
time according to (31), It Is not possible to rigorously calculate 
the corresponding eigen values. An approximate* calculation must 
be used. In order to avoid special rod dimensions. It Is appro- 
priate to make the system (32) and the co’’ res ponding boundary 
conditions (33) dimensionless. In addition to the position coor- 
dinate 



we will Introduce the nondlmenslonal variables 


the stiffness ratios 


and the eigenvalue 



T mr, 


Eh 


A 


Ml 


“ Eh 


(35) 

(36) 

(37) 
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If we use the new coordinates and parameters, we find that 

V- _ 0 - i j r + (, - i) r]" ... 1 ^ ^ 

J r- + (1 - »t/" - -j- 1 1 /' - - - j) w]' - • , 

with the boundary conditions 

t'mmo, U" + ~ m 0 , f •• O , ffMiOnndl* (39) 

This dimensionless elrenvalue problem can then be associated 
with the variation problem, so that the connection with one of 
the direct methods of variational calculus has been establ 1 shed . 



In order to approximately calculate the olpenvalues using the 
Rltz method, the variational problem (40) 1s stated In algebraic 
terms. For this purpose, we w’ll use a single term trial solution 
for the calculation discussed here 

r(n) - 

We will use the eigen functions 

« r, « »in .Tff (42) 

as position functions, which ar-e the elcen ^unctions of the tipping 
support for which the nonlinear parts In the torsion q (14) can be 
Ignored [4], These satisfy all of the geometric boundary conditions 
In (39) as required. 



Figure 3: Critical tipping moment of a rectangular 

profile support. 

After carrying out the variations given In (40), 


cross-section 
we obtain a 
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homogeneous equation system For ^ free ooefflclents a^, In 
(^1). If the associated determinant. Is set equal to 0, then v<1'\ 
( 36 ) and (37) we find the corresponding critical tipping momenv 
with the lowest order to br 


If, as In [^4], 


A/ 

we had 


jl 

, «/» •cyr ‘ 
Ut~ iKJ, 
Ignored the rod 


Inclination during the 


(^ 3 ) 


differentiation, we would have obtained 


On the other hand, 
have obtained 


M 


Kt 



from rigorous 




I 


*7, 

~^J, 

nonlinear equations we would 





(^5) 


[9]. ‘jr EJ 2 KJj and GJ,p<<EJj^, all of the formulas (^<3) to 

(^5) lead to the Prandtl tipping moment 

(46) 

We will qualify our result using the example of a rectangular 
steel strip (Height h, width b, G * O. 385 E). The nond Imenslonal 
tipping moments * 1» 2, 3) are plotted as a function of 

the width-height ratio b/h In the range 0 < b/h ^ 1/2, feee Figure 
3). It Is known that the critical moment can Increase considerably 
with Increasing width. Comparison of the result shows the useful- /i(,q 
ness of the approximation (43) and (44). It seems that (44) is 
even more suitable In practice, because In this way we have a 
certain safety margin. By using a multiple-term trial solution 
(4l) It would be possible to Improve the result (43). From the 
theoretical point of view, (43) la better than (44). 


5. SUMMARY AND OWRVIEW 

Tipping problems of elastic rods are of theoretical and prac- 
tical interest. Usually, one has to deal with complicated boundary 
value problems. This means that It Is advantageous to derive the 
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differential equations and boundary conditions from a variational 
principle. 

One of the main assumptions In this process up to the present 
was a very small deviation between the undeformed and deformed 
rod axis, so that differentiation with respect to the position 
coordinates could be considered to be Identical. Here this re- 
striction Is disregarded. If the central line does not experience 
any longitudinal extension, then froTi the analogy between the motions 
of a gyroscope and the elastic behavior of a rod, one can derive 
nonlinear relationships for the curvatures and for the torsion. 

These compliment previously known relationships because non- 
llnearltles occur In the torsion which are of the same order of 
magnitude as In the other characteristics of the rod axis. 

Using the method of small oscillations of small expanded 
boundary value problem, we derive the corresponding stability 
equations. Using the simple example of a bendlrg-tippl ng problem 
we demonstrate the Improved perfcrmance of the suggested equations. 
Among other things, we find that the constant direction moment 
and < he following moments are not conservative [10], a fact 
whlr.h Is known from twisted rods. 

If It Is desired to extend our results to the compressed 
tipping rod, then a nonlinear longitudinal extension will be 
added to the previously mentioned nonlinear quantities. We may 
assume that the expressions (14) and (15) for p^ ^ and q will 
remain unchanged and that the expression 

(«■• 4 - 4 k]r't) 

will be used for the extension of the rod longitudinal axis. Just 
like In [4]. 
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